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ABSTRACT 

This paper presents the design and evaluation of 
three regulators for an industrial evaporation 
station: a conventional PI Controller, a Fuzzy Logic 
Controller and a Generalized Predictive Controller 
(GPC). Linearized models are derived to design the 
predictive controller. The process diagram, the 
instrumentation, the nonlinear model, the 
controllers design and the implementation are 
described. Performance of the three regulators are 
compared by simulation. 



1. INTRODUCTION 

One of the main stages in the sugar 
production, takes place in the evaporation station, 
where the sugar content of a juice rich in saccarose 
is increased by boiling. The resulting syrup is used 
to obtain sugar crystals in a set of vacuum pans. 
For the study we considere a specific beet sugar 
factory placed in Valladolid (Spain), which 
evaporation section is composed of five 
evaporators working as a four-effect system. 

One of the most important aim of the 
control system of an evaporation station is to 
maintain the sugar concentration at the output of 
the station at prescribed values despite the 
disturbance effects. It is also very important to keep 
the level of juice in each evaporator between two 
security limits. 

The complexity of the system dynamic and 
the presence of disturbances justify the use of 
advanced controllers. 

The paper describes first the plant and the 
instrumentation. A derivation of a nonlinear model 



of the station follows. Next, the plant is analized by 
simulation. The design of two regulators, based on 
fuzzy logic and generalized predictive control is 
then presented. The performance of both regulators 
is compared by simulation of the nonlinear 
dynamical model with the performance of a 
conventional PI controller. Conclusions are 
presented in the last section. 



? PLANT DESCRIPTION 

Fig. 1 shows a two chambers evaporator. 
The heating chamber surrounds a set of vertical 
tubes that contain boiling juice. A flow of steam 
enters to this chamber and transfers heat to the 
juice providing the energy needed for boiling. 




Fig. 1 : A two chambers evaporator 



The steam condenses around the tubes 
and leaves the evaporator as condensate. The 
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interior of the tubes, plus the evaporator upper and 
bottom spaces, form the juice chamber. A sugar 
solution of low concentration (juice) flows 
continuously into the base of the evaporator and 
starts boiling. Consequently, a solution of higher 
concentration is obtained at the output. The steam 
produced from the water evaporation reaches the 
upper space and leaves the juice chamber by a 
pipe at the top. 
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Fig. 2: Process diagram 

Fig. 2 shows the process diagram. It 
consists of a set of five evaporators interconnected 
through pipelines and valves. The steam generated 
in one evaporator is used to provide energy to the 
next heating chamber, while the juice flows from 
one evaporator to another increases the sugar 
concentration. In this multiple-effect arrangement 
only the first evaporator is fed with fresh steam and 
juice. In the last evaporator, the evaporated steam 
scurries from the juice chamber to the condensers 
and then to atmosphere. Since the pressure in the 
last heating chamber is very low, a set of 
condensers is required to decrease the vapour 
pressure in the juice chamber too. In this way the 
temperature at which the juice starts boiled drops to 
a lower value and evaporation is possible. 




Fig. 3: Instrumentation and control structure 



Fig. 3 shows the instrumentation and 
control strategy. The control strategy includes four 
level indicator controller, two preassure indicator 



controller and a cascade concentration-flow control 
system. The instrumentation includes also 
measurement of the vacuum pans demand. 



3. NONLINEAR MODELLING 

The nonlinear mathematical model of the 
plant is obtained by applying simple physical laws. 
First we consider the modelling of a single 
evaporator. Later, we link all the models together. 

Balances of mass, sugar and energy of the 
juice chamber of a evaporator lead to the following 
equations at time t: 



dm/di = q t - q - q e 
d(mc)/di = qj c t - q c 



(1) 



(2) 



d(mh)/dt = qjhj + Q-qh'q e H (3) 

where q- t and q are the input and output mass flows 
of juice respectively, c,- and c are juice sugar 
content of both flows and q e is the steam flow 
obtained by boiling. The juice enthalpy h is a known 
function of c and temperature T. Steam enthalpy H 
depends on pressure p and temperature T. 

As hypothesis we assume that the steam in 
the vapour space of the juice chamber is in 
equilibrium with the boiling juice. Then, the 
pressure in this chamber is determined by the 
equilibrium relation linking T t c and p. 

The heat from the steam to the juice is 

Q = K(T S -T) (4) 

where the heat transmission coefficient K is a 
function of the juice level in the tubes / and c. T s is 
the saturation temperature. 

The time constants associated with mass 
and energy balances in the heating chamber are 
neglectable compare to those in the juice. So, we 
can assume steady-state behavior, i.e.. 

Q = F s (H-h) (5) 

log(P s ) = 2147 / <T s +273.2) + 5.76 (6) 

F s is the steam flow that enters the heating 
chamber. H-h is the latent heat of steam which is a 
function of the saturation temperature T s (°C), 



649 



related to its pressure P s (bars) by Antoines law 
(6). 

If we also assume steady-state in the 
vapour space of the juice chamber, then 

Qe = K v ^(P* * Po*> & 

K v being an installation dependent coefficient and 
p 0 the atmosphere pressure. 

Finally. q\ and q are functions of valve 
openings as well as of pressure and level of juice 
(measured from the valve position): 

q =K, sqrt((p + dj g I - Po) dj /(1 '« 1 2 + K j1 )> < 8 > 

q/ =K;2 sqrt((Pt -djgl-p) dj I{Vu# + K ;3 )) (9) 

Ui and u 2 are valve openings and the Kjs are 
valve dependent coefficients; p| is a contour 
pressure. The level / can be obtained from the total 
juice mass m, level juice t Q% mass juice m Q (at the 
bottom of the tubes) and the evaporators cross 
section surface S as 

f=l 0 +(m-m 0 )HSdj) (10) 
The juice density dj can be calculated from 

T and c. 

The model of the whole evaporation station 
is obtained by repeating all the equations five 
times, and equating outputs flows and pressures of 
an evaporator with the corresponding input of the 
next one. A library of physical properties of steam 
and sugar juice is used to obtain the juice and 
steam enthalpies, densities, etc. Other data as 
valve coefficients. installation dependent 
coefficients and so on are taken from the technical 
documentation of the factory. 



A. PLANT ANALYSIS 
4.1 Control problems 

In this process there are two main control 
aims to consider 

• Maintaining steady constant concentration of 
sugar at the output, specially to reduce large 
and repetitive disturbance effects. 

• Keeping the juice level in each evaporator at 



prescribed values for security. 

We use PI control to regulate levels of fluid 
and pressures (see Fig. 3). A cascade PI controller 
regulates the juice sugar content and the juice 
flowrate. Note that variations in the juice flowrate 
change the resident time of the juice within the 
evaporator and consequently the concentration. 

To maintain the level of juice, the control 
strategy varies the valves opening placed in the 
juice line at the input of each evaporator as it is 
shown in Fig. 2. Pi control is sufficient to achieve 
this objective. 

The most severe disturbance is the varying 
steam demand from the vacuum pan station which 
influences the operation of the system and the final 
output quality. The concentration is not more as 
constant as desired. 



Fig. 4 shows a typical disturbance and its 
effect on the flowrate of juice at the input of the 
tank. 
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Fig. 4: Steam demand and output juice 



concentration 

In order to design a concentration regulator 
providing good control performance we need a 
control system that results in a good disturbance 
rejection and robustness against the non-linearities 
present in the system. Fuzzy logic and GPC based 
regulators seem to be good choices for our 
purposes as we will see in the next sections. 

4.2 Open-loop responses 

Fig. 5 shows open loop response of the 
juice concentration, when applying a negative step 
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in the flow set-point. 

It could be seen, clearly, the existence of a 
delay in the process and the typical response of a 
non minimal phase system. 
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Fig. 5: Control signal and open loop response 

4.3 Linear modelling 

By using the commercial software ACSL is 
possible to obtain from the nonlinear model, 
simplified linear models valid nearly specific 
operation points. We work out two function 
transfers in order to describe the relationship 
between output concentration c, flowrate-reference 
q$* and steam demmand p^. MATLAB was used to 
reduce the order and discretize the models with 
sampling time of 0.06 hour. The transfer functions 
are the following: 

Output flow rate set point 

^5 <f 1 (' 0. 127 + 0. 13 CT 1 - 0.048 Cf 2 ) 
q 5 * 1 - 1.8 cf 1 + 1.055 <f 2 - 0.184 <f 3 



(11) 



(12) 



Steam demand 

tf 1 (3. 74 + 3.93 cf 1 - 0.19 cf 2 ) 
p d 1 - 1.8 cf 1 + 1055 cf 2 -0.184 cf 3 

5. REGULATORS DES1GM 
5.1 Fuzzy logic control 



Fuzzy logic provides an alternative solution 
to nonlinear control. A typical fuzzy expert 
controller comprises three elements: membership 
functions to fuzzify the physical inputs, an inference 



engine with a decision rule base and a defuzzifier 
that converts fuzzy control decisions into physical 
nonfuzzy control signals. With these components, a 
complex nonlinear controller can be built. 



Knowledge 
base 



Inference 
making 



DefuzzHter 



PROCESS 



Fig. 6: Structure of the fuzzy logic regulator 

Fig. 6 shows the block diagram of the fuzzy 
controller used in this study. The inputs to the fuzzy 
controller are error e(t) and 



de(t)~e(t) - e(t-1) t 



(13) 



the difference between the present and the 
previous error. The output du(t) is the incremental 
change of the control signal given by 



u(t) = u(t-1) + du(t) 
GE, GR and GU are parameters to be tuned. 



(14) 



Table I presents the incremental control 
rules, in form of a fuzzy control matrix. This fuzzy 
control matrix, proposed by Tang (1987), involved 
seven linguistic terms to characterize each 
variable: NB (negative big), NM (negative 
medium), NS (negative small), ZE (zero), PS 
(positive small), PM (positive medium) and PB 
(positive big). 

The control matrix is based on the following 
principles: 

• If the output has the desired value and the error 
change is zero, then the output of the controller 
remains constant {du(t)-O). 

• If the output y(t) diverges, the action depends 
on the signum and the value of error e(t) and 
the error change deft). If the conditions are 
such that the error can be corrected quickly by 
itself, the controller output remains constant or 
almost constant. Otherwise, the controller 
output is changed to achieved satisfactory 
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results. 
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Fig. 7: Membership functions 

As Fig. 7 shows, in this study we use 
triangular membership functions. An standard 
center of area method was used for defuzzification. 



5.2 Generalized predictive control 

A generalized predictive controller is based on a 
model of the system decribed by 

A y(t) = cf 1 B u(t) ♦ C v(t) * D/A w(t) (15) 

where y(t) is the output, u(t) the input, v(t) a low 
frequency disturbance, w(t) a white noise, A, B t C 
and D polynomials on the delay operator q"\ and 
A= 1 -q'*. 

The predicted output ytft+j) over a range 
up to the prediction horizon can be obtained as a 
function of the polynomials A t B, C and D, the past 

outputs ]y(t).y(t-D 1. ™ ure ** 

[r(t+V,r(t+2)..l past controls \u(U1),u(t-2h \ aruJ a 
vector of potencial future controls [u(t) t u(t+1), . 

The GPC algorithm calculates the vector of 
future controls by the minimization of the cost 



function: 

Jcr,- = X c2 C + °M' + 7 " »-"(' + > " 2)] 

(16) 

where n1 is the "costing horizon", n2 is the 
"prediction horizon", nu is the "control horizon" after 
which the control increments are supposed to be 
zero, a is a factor that weights the control efforts. 
The error eft) is given by 



e(t+j)=r(t+j)-y p (t+j) 



(17) 



We can get feedforward actions to 
compensate for the disturbances v(t) by including 
the C polynomial in the model. The noise affecting 
the output of the system is also consider in the 
control design through the polynomial D of the 
model which can be modified to improve the 
control performance. 



ff cumin ATinN EXPERIMENTS 

6.1 Basis for the evaluation 

A set of simulation experiments with 
differents controllers are presented in this section. 
The output of the system is the concentration (C5) 
and the control signal is the flowrate setpoint (Q 5 *>. 
The disturbances acting on the system are the 
steam demand from the vacuum pans (p<j) and the 
pressure set point in the juice chamber of the fifth 
evaporator (P5O (see Fig. 2). 

6.2 PI control 

This experiment was carried out by using a 
continuous PI controller to maintain the 
concentration close to the reference signal. The 
controller was tuned by a trial and error procedure 
that gives following values for the proportional gain 
K p and the integral time Tf 



K p = -3 



T/ = 0.02 hour 



Fig. 8 shows the set-point C5*. the output 
c 5 and the control signal 95*. As we can see the 
variance of the concentration is quite high mainly 
because of the disturbances acting on the system. 
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Fig. 8: Plant response and control signal with 
PI control 

6.3 Fuzzy logic control 

The regulator used in this experiment 
corresponds to the controller described in 5.1. The 
controller is characterized by the following 
parameters: 



GE = -1 GR = -0.05 GU = 30 
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Fig. 9: Plant response and control signal with 



fuzzy control 

Fig. 9 shows the set-point c 5 *. 
concentration c 5 and the control signal q$* 
obtained with fuzzy control. We appreciate that 
concentration response improves considerable, 
compared with Fig. 8. 

6.4 GPC control 

We implemented first the control algorithm 
without considering the feedforward polynomial 



(C=0). The design parameters in this case were: 

n1 = 2 n2 = 10 nu = 2 a = 0.03 

The results by considering constant 
reference are given in Fig. 10. The process 
behaviour with GPC control is also better than with 
Pi control. 

In order to reduce the variance of the 
output C5 we consider in the model the C 
polynomial corresponding to the steam demand. A 
feedforward component is therefore included in the 
control law. Fig. 11 shows that the output variance 
has been drastically reduced. 

The influence of the weighting factor on the 
control efforts can be observed by comparing the 
last result with Fig. 12, obtained with a=0.01. In this 
case the control efforts are greater and the output 
variance is smaller. 
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Fig. 10: GPC without feedforward action 
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Fig. 11 : GPC with feedforward action 



653 



rrrr 

*>^> 


Concentration c5 


i 


t- 








time (min) 


. „ . _j 




t 

r ~~ r '" " ' 

! 


! ■- t 




i 
i 

time (min) 





Fig. 12: Effect of the weight factor 



6.5 Comparison 

Table II presents a quantitative comparison 
of the performance of the three regulators, based 
on experimental results presented on Fig. 8, 9 and 
11. Mean value and root mean square of two 
important variables, error e(t) and incremental 
change du(t), are considered as comparison index. 

Table II 



Indexes of 
peffofmence 


PI 

Control 


Fuzzy 
Control 


OPC 
Control 




0.0015 


0.0025 


0.0045 


rms-e 


0.3179 


0.2362 


0.2628 


meen-du 


0.0009 


0.0012 


0.0066 


rrnS'du 


0.5144 


2.7287 


2.0438 



for purposes of implantation in industrial plants. 
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7. CONCLUSION 

The work carried out allowed to evaluate 
two advanced control techniques applied to a 
process of complex dynamic behavior as the 
evaporation station. 

The results confirm the better performance 
of fuzzy and predictive controllers over PI control. 
The juice concentration shows much less variance, 
in spite of the disturbance that it affects. 

The fuzzy controller doesn't require a 
mathematical model and the computational 
requerimerrts are less in comparison with GPC 
controllers. These differences result quite attractive 
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